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Constructing qubits which are suitable for quantum computation remains a notable challenge. 

Here, we propose a superconducting flux qubit in a dc SQUID structure, formed by a conventional 
insulator Josephson junction and a topological nanowire Josephson junction with Majorana bound 
states. The zero energy Majorana bound states transport 4ir period Josephson currents in the 
nanowire junction. The interplay between this 47r period Josephson effect and the convectional 2tt 
period Josephson effect in the insulator junction induces a double-well potential energy landscape 
in the SQUID. As a result, the two lowest energy levels of the SQUID are isolated from other levels. 

These two levels show contradicting circulating supercurrents, thus can be used as a flux qubit. We 
reveal that this flux qubit has the merits of stability to external noises, tolerance to the deviation of 
system parameters, and scalability to large numbers. Furthermore, we demonstrate how to couple 
this flux qubit with the Majorana qubit by tuning the junction parameters, and how to use this 
coupling to manipulate the Majorana qubit. 

PACS numbers: 74.90.+n, 74.50.+r, 03.67.-a 


I. INTRODUCTION 

Quantum computation, which involves processing in¬ 
formation on quantum variables, is valuable in directly 
simulating quantum mechanics or efficiently solving clas¬ 
sical problems^. The realization of the quantum compu¬ 
tation calls for the coherent storage and manipulation of 
the quantum information^. A typical scheme for coherent 
quantum information processing consists of two energy 
level systems, the so-called quantum bits (qubits)^. In a 
qubit, the two levels serve as the basis states, then the 
quantum superposition of these basis states can store the 
quantum information. A number of qubits have been pro¬ 
posed in both natural and artificial two-level systems^ - —. 
In particular, the artificial qubits seem to be very promis¬ 
ing for realistic application. In fact, more than one hun¬ 
dred artificial qubits have been integrated in the recent 
commercial quantum computers^. 

One attractive candidate among the artificial qubits 
is the superconducting flux qubitiii - — , which is com¬ 
posed of a superconducting loop interrupted by one or 
more Josephson junctions. Flux qubits have several out¬ 
standing advantage s 11 ! 12 . The first one is the scalability: 
flux qubits are based on lithography technology which 
allows a quantum computer to integrate a large number 
of qubits. Second, flux qubits are resistent to the setback 
charge fluctuation, which is the major noise in the elec¬ 
tronic systems. The last and most important advantage 
is the flexibility: the property of the the flux qubit is de¬ 
termined by the design of the Josephson junction system, 
in which the parameters can be tuned continuously. Two 
designs for flux qubits have been studied up to nov&^: 
(i) an rf SQUID with one Josephson junction; and (ii) 
a superconducting ring with three Josephson junctions. 
The one-junction design demands a high self-inductance 
of L - 200pH, which can be realized only in a large 
device. For this reason, the qubit is usually very sus¬ 
ceptible to external noises^. The three-junction qubit 


is constructed in a superconducting loop of micrometer 
size with a very small self-inductance L ~ 5pH and low 
critical currents I c ~ 200nA. Therefore, this qubit has 
weak coupling with external fluctuation fields and was 
adopted more widely^. In this three-j unction design, 
however, the parameters of the junctions are severely 
confined. Two of the three junctions must have exactly 
the same critical current while the third junction has a 
slightly smaller critical current. This strict constraint 
makes the three-j unction design difficult to be realized 
in experiment s 11 1 12 . The deficiencies of the one-junction 
and the three-j unction flux qubits urge the development 
of a third type of flux qubits with both low inductance 
and weak restriction on the junction parameters. This 
might be accomplished by exploiting exotic Josephson 



FIG. 1. (Color online). Schematic setup of a flux qubit in 
a dc SQUID, formed by a conventional insulator Josephson 
junction and a topological nanowire Josephson junction. The 
wire hosts four Majorana bound states 71,2,3,4. The two ar¬ 
rows indicate the two contradicting circulating directions of 
the supercurrents, which are the basis states of the qubit. 
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junctions. 

Recently, the topological nanowire Josephson junctions 
have drawn much attention because the zero energy Ma- 
jorana bound states (MBSs) appear in the ends of the 
wires!! - —. The MBSs are potentially useful for fault- 
tolerant quantum computation: two MBSs can form one 
Majorana qubit, and the braiding of the MBSs can ro¬ 
tate the Majorana qubit!! - —. In Topological junctions, 
MBSs show several unique features which can be detected 
in experiments!! - ! 5 . One important feature of MBSs is 
that they transport 47r period super currents, known as 
the 47 t period fractional Josephson effec t 15 i 22 i 26 . This 
novel Josephson effect is regarded as a compelling evi¬ 
dence for the existence of MBS o 18 i 22 i 26 . Besides, it can 
be used to couple the Majorana qubits and the super¬ 
conducting qubits coherently!!. This coherent coupling 
can construct quantum gates for Majorana qubits, which 
is important for the quantum computation!! - —. The 47 t 
period Josephson effect has been experimentally reported 
in topological nanowire Josephson junctions, where indi¬ 
rect ac signals were observe d 32 ! 33 . The direct experimen¬ 
tal observation of the fractional Josephson effect is hin¬ 
dered by several factors. The first one is the existence of 
a competing 2 tt period supercurrent from other Andreev 
bound states!! - —, which might be even larger than the 
47r period component. In addition, the small coupling 
between distant MBSs will induce a parity flipping in 
the topological junction, which diminishes the 47r period 
signal in the dc current!!. 

In this work, we propose a two-junction flux qubit us¬ 
ing the 47 t period Josephson effect. For this purpose, 
we adopt a setup as sketched in Fig. 1 . A conven¬ 
tional superconductor-insulator-superconductor Joseph¬ 
son junction is connected with a topological nanowire 
Josephson junction, forming a small SQUID with van¬ 
ishing self-inductance. In the middle of the nanowire 
junction, a tunneling barrier is created by an external 
voltage gate. Four MBSs exist in the wire: 72 and 73 
locate near the tunneling barrier, 71 and 74 locate at 
the ends of the wire. The tunneling between 71 and 72 
contributes a 47r period component to the Josephson en¬ 
ergy of nanowire junction. This 47r period Josephson 
effect induces a double-well potential in the system, un¬ 
der which two lowest energy levels are isolated from other 
high energy levels. These two isolated levels present con¬ 
tradicting circulating super currents, thus making up a 
flux qubit. We study the the system with both analyti¬ 
cal and numerical methods, and find that the flux qubit 
exists in a wide junction parameter range. Therefore, 
our proposal provides a third type of flux qubit with low 
inductance and large variability in junction parameters, 
which bring in the advantage of resistance to the external 
noise and robustness to system parameter modulations. 
Moreover, we show that this flux qubit can be coupled 
with the Majorana qubit when the junction parameters 
are well tuned. This coupling is useful for manipulating 
the Majorana qubit, which is important for the topolog¬ 
ical quantum computation. 


This work is organized as follows. We present the 
model in section II, and simplify the the system to an 
equivalent one particle Hamiltonian. In section III we 
study the potential energy of the system, and shows that 
double-well structure exists. In section IV we discuss the 
quantum tunneling within the double-well potential and 
analyze the parameter regime for a flux qubit. In section 
V, we show the numerical simulations on the eigenstates 
of the system, and discuss the energy structures of the 
flux qubit. In section VI, we discuss the coupling between 
the flux qubit and the Majorana qubit. Finally we give 
discussions and a conclusion in section VII and VIII. 


II. MODEL 

The setup shown in Fig. 1 is built by a conventional 
insulator Josephson junction and a topological Josephson 
junction. The topological junction is constructed with a 
spin-orbit coupling nanowire between two conventional 
s-wave superconductors. The nanowire enters the topo¬ 
logical superconducting phase with appropriate Zeeman 
energy and chemical potential. In the middle of the wire, 
a voltage gate is applied to produce a tunneling barrier. 
Four MBSs exist in the wire, two of them locate at the 
ends of the wire and the other two locate near the tunnel¬ 
ing barrier. The low energy Hamiltonian of the nanowire 
junction is determined by the the superconducting phase 
difference across the junction and the quantum states of 
the MBSs!!*!!*!!*!!, 

Q 

Um =-E c idj m - J m cos 0 m + iE m cos -^7273 

+«<5l7i72 + *<5r7374, (1) 

where 72 and 73 are the two MBSs near the tunneling 
barrier, 71 and 74 are the MBSs at the ends of the wire, 
0 m is the phase difference, E c \ = e 2 j 2 C\ is the charging 
energy from the capacitance C\ of the nanowire junction, 
Sl and Sr represents the small wave-function overlaps 
between the two MBSs at the left and the right sides 
of the tunneling barrier, respectively. J m and Em ac¬ 
count for the 27 r and 47r period Josephson energy from 
the quasi-particle channels and the MBS channel. The 
Jm might be larger than E m , since we have large num¬ 
ber of quasi-particle channels and only one MBS channel. 
This produces one of the obstructions for the successful 
detection of the 47 t period Josephson effect. We note 
that the 47r period Josephson energy will induce an 2 tt 
period dc Josephson current after adding the MBSs cou¬ 
pling 5 l,r- In this sense, a 47 t period dc supercurrent 
does not exist. The 47r period Josephson effect actually 
refers to the Josephson energy shown in Eq. ( 1 ). 

The conventional insulator Josephson junction has 
been extensively investigated, with a low energy Hamil¬ 
tonian given as!!, 

Hj = -E^Oqj — J cos #j, 


( 2 ) 
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where Oj is the phase difference, E c 2 = e 2 /2C2 is the 
charging energy, and J is the Josephson energy. A mag¬ 
netic flux is applied through the SQUID, which induces 
a inductive energyii, 

U L = - 0 m - 2tt$/$ 0 ) 2 , (3) 

where L is the self-inductance of the SQUID, 4> is the ap¬ 
plied flux through the loop measured in the unit of mag¬ 
netic flux quanta <f>o = h/2e. In our setup, we hope to 
minimize the environmental magnetic noises. Therefore, 
we adopt a small SQUID with a vanishing self-inductance 
L « 0. In this scenario, the inductive energy should al¬ 
ways be minimized by a constraint^ between 0 m and #j, 

e = e rn = e J - 2 tt $/$ 0 . ( 4 ) 

Under this constraint, only one phase difference 0 is 
enough to describe the system, drastically simplifying the 
problem. Adding together, we arrive at a total Hamilto¬ 
nian for the system, 

T-L =—E c 8q — J cos (0 + 27t<I>/4>o) — J m cos# 

Q 

+iE m cos -7273 + ^l7i72 + ^r7374, (5) 

where E c = E c \ + E c 2 is the total charging energy. 

The four MBSs in the nanowire form two Major ana 
qubits which reflect the parity of the superconducting 
ground state of the nanowire. The topological parity 
states is clearly presented by defining complex fermionic 
operators with the Majorana operators, 

fi = (74 - *71)/2, A = (72 - *73)/2. (6) 

These two fermionic operators are equivalent to the four 
MBS. Their eigenstates |0), fl |0), /JjO), and /1/2 |0) cor¬ 
respond to the superconducting ground states with dif¬ 
ferent fermionic parities. The total parity of the system 
should be conserved at low temperature, thus we can de¬ 
compose the Hilbert space of the Hamiltonian ([5]) into 
two disconnected sub-spaces with even or odd total par¬ 
ity. We take the even total parity without losing gen¬ 
erality, then the subspace is expanded by |0)i|0)2 and 
/t/t|0>. It is instructive to reformulate the parity states 
into pseudo-spin states ||) = |0)i|0)2 and |t) = flfl |0). 
In this representation, the Majorana operators can be 
transformed into pseudo-spin operators 27273 — cr z and 
27172 = —^7374 = cr Xl and the Majorana qubits are trans¬ 
formed into spin qubits. Now the Hamiltonian is rewrit¬ 
ten as, 


and a magneto-static potential cos | + S m a x . We 

can write down the Hamiltonian explicitly in a 2 x 2 form 
using the spin states in z direction as the basis, 


n = 


€ + u + ( 0 ) 2 S m \ 

^ + u-(e)) , 


(8) 


where Pq = ihde is the momentum operator for 0 , 
£/*(#) = — J cos (0 + 27 r 4 >/ 4 >o) — Jm cos 0 =b Ern cos | cor¬ 
responds to the potentials for the up- and down-spins 
respectively, where U + ( 6 ) = U ~(0 + 2 tt). The off- 
diagonal term S m accounts for the spin flipping strength 
from the wave-function overlap between MBSs, which is 
small compared with the Josephson energy J for a long 
nanowire. The same effective Hamiltonian can be ob¬ 
tained for odd total parity sub-space by simply taking 
Sm = $L + Sr. 


III. DOUBLE WELL POTENTIAL 

As pointed out in previous studie s 11 ! 12 ! 37 , a double-well 
potential is the key ingredient for isolating two energy 
levels in the one particle Hamiltonian established in Eq. 
(8). To obtain a double-well potential, we consider half 
a quanta of applied magnetic flux through the SQUID. 

= 4>o/2, then the potential U ± (#) can be written as, 

Q 

U ± ( 0 ) = E m (a cos 0 =b cos -), ( 9 ) 

where a = (J — Jm)/E m - We find that the 27 r period 
component J m in the topological junction is canceled by 
J in the conventional junction, through destructive in¬ 
terference induced by the magnetic flux 4 >. When J is 
tuned close to J m , the Air period component E m will 
make significant contribution to the total Josephson en¬ 
ergy, even though it is smaller than J m . We show the 
potential U ± { 6 ) for a = 0.2, 0 . 5 , and 1.5 in Fig. 2 a, 
2 b and 2 c, respectively. The potential for small a has 
no double-well structure as shown in Fig. 2 a . In Fig. 
2b, we immediately find the double-well potential struc¬ 
ture, which is required for a flux qubit. Figure 2 c in¬ 
dicate that the double-well potential still exists when a 
becomes larger, however, the potential barrier within one 
double-well increases and approaches the tunneling bar¬ 
rier between different double-wells. We hope to clarify 
the parameter regime for the double-well potential struc¬ 
ture. Therefore, we study the classical stability points 
which are given by the zeros in the first derivative of the 
potentials, 


H =—E c dl — J cos (0 + 27t4>/4>o) — J m cos 0 

Q 

+cr z E m cos - + 

Srn &x-> ( 7 ) 

where S m = $l — Sr. The Hamiltonian is now equiva¬ 
lent to a spin one half particle using 6 as the coordinate. 
The particle has a mass of m = h 2 / 2 E C , moving under 
electro-static potential — J cos (6 + 27 r < h/ 4 >o) — J m cos# 


517 * ( 0 ) 
30 
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0 1 


2 E^ sin - [a cos - =b - = 0. 


( 10 ) 


We immediately obtain two sets of solutions for this equa¬ 
tion when a > 0: 0 = 2 nn with an integer n; and 
0 = ±2 arccos (d=l/ 4 a) + 4n7r only for a > 1 / 4 . Now the 
parameter space is divided into two regimes by a critical 
value a c = 1 / 4 . In the regime of a > <a c , we have two 
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FIG. 2. (Color online). Potential functions U + (0) (solid line) and U~ (0) (dashed line) as a function of superconductor phase 
difference 0 , with (a) a = 0.2, (b) a = 0.5, and (c) a = 1.5. (d) Schematic illustration of three nearest neighbour quantum 
tunneling process: (i) spin-flipping, (ii) out of the double-well, and (iii) within the double-well. 


minimum points and two maximum points, which is the 
double-well structure as shown in Fig. 2 b and 2 c. In 
the regime of a < <a c , there are only one minimum point 
and one maximum point as shown in Fig. 2 a, then the 
double-well structure disappears. For a < 0 , we would 
have the potential up-side-down, and double-well struc¬ 
ture does not exist. 

In our proposal, we exploit the destructive interference 
in the SQUID to reduce the 2 i r period component in the 
Josephson energy, and highlight the 47T period compo¬ 
nent. We note that the Josephson energy of the conven¬ 
tional junction plays a two-fold role in this setup. First, 
it cancels the 2 tt period components in the Majorana 
junction. Second, it interplays with the 47r period part 
to create the desired double-well potential energy land¬ 
scape. 


IV. QUANTUM TUNNELING 

We want to construct a flux qubit, thus take the regime 
of a > a c where the double-well structure exists in 
the potential U ± ( 0 ). We consider the local states at 
the minimums of the potential, and analyze the energy 
spectrum due to the quantum tunneling between adja¬ 
cent local minimums, within the standard tight binding 
formalism^!. In usual double-well potential systems, the 
two bonding and anti-bonding states should be isolated 
in the double well, where the probability of the quan¬ 
tum tunneling out of a double-well is negligible. In this 
system, however, one subtlety must be addressed: the 
particle has a spin, while the spin-up and the spin-down 
particle have different potentials. As a result, we have 
to consider the quantum tunneling from the spin-flipping 
process induced by In the following, we study all 
three types of nearest neighbor quantum tunneling as 
sketched in Fig. 2 d: (i) the spin-flipping quantum tun¬ 
neling between the nearby local potential minimums for 


different spins, (ii) the tunneling out of a double-well 
in one spin branch, and (iii) the tunneling within one 
double-well. We will calculate the three hybridization en¬ 
ergies from each quantum tunneling processes, and find 
the parameter regimes in which two of them can be ne¬ 
glected. 

We start from the spin-flipping quantum tunneling be¬ 
tween the nearby local states of different spin branches. 
The hybridization energy is determined by the overlap 
of the two local wave functions, multiplied by the spin¬ 
flipping energy 2 6 m , 

Ex « 2S m x(8x,02) = 2 5 m j d9^*{9,9^(9,9 2 ), ( 11 ) 


where 0 \ = 2 arccos ( 1 / 4 a) is a local minimum for U~ ( 0 ), 
and 62 = 27 r — 0\ is the local minimum for U + ( 6 ) close to 
01; -0(0, #1,2) are the lowest energy local wave-functions 
at the potential minimum 0i ? 2, which can be calculated 
with a harmonic oscillator approximation. We define the 
harmonic potential V(6) = \K (0 — $1,2) 2 + Vo, where 
Vo = ~E rn (a + 1 /8a) corresponds to the minimum of 
F r± ( 0 ), and K is given by, 


K = 


d 2 U~( 0 ) 

QQ2 l»=*l 


Em (o? 



( 12 ) 


The local wave function is obtained by solving the 
Shrodinger equations under the harmonic potential, 


^ 6i) = ( 2 ' (13) 

with i = 1 , 2 . Then the wave-function overlap is given 

by, 


X{9i,9 2 )=e 02)2 . ( 14 ) 

Plugging this wave function overlap back to Eq. ( 11 ), 
we can obtain the hybridization energy from the spin¬ 
flipping tunneling. 
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FIG. 3. (Color online). Hybridization energies E\ (solid line), 
E 2 (dashed line), and E 3 (dotted line) as a function of a, with 
E c — 0.05 Em and 5 m — 0.0001F m . 


Second, we consider the hybridization energy from the 
quantum tunneling out of a double-well potential. We use 
the WKB method to solve the problem, which was widely 
used for treating the quantum tunneling in Josephson 
junction systems^!. The resulted hybridization energy is 

ujh ( 1 r^ir-Oi \ 

(15) 

where uo = \/2KE c /h is the frequency of ground state 
within the harmonic oscillator approximation. The in¬ 
tegration on the exponential function represents the ac¬ 
tion of the particle moving from one minimum potential 
point to another. With the same formalism, we can also 
obtain the hybridization energy from the quantum tun¬ 
neling within one double-well, 

E 3 « ^exp (^-y=F /; dffy/U-(e)-V^j . (16) 

We notice that the formulas for the two hybridization 
energies E 2 and E% are identical, except different actions 
from the different integration zones. The value E 2 /E 3 
is depending on the potential landscape, which is con¬ 
trolled by a. We expect that E 2 /E 3 should approaches 
one for large a and approaches zero for small a. E\ is 
proportional to S m which is an exponential function of 
the nanowire length. Therefore, E\ will also be very dif¬ 
ferent for long and short nanowires. Two of these three 
energies should be vanishingly small to construct a flux 
qubit, i.e. E% E \ ? 2 in this system. The result of ^ 1 , 2,3 
is shown in Fig. 3, as a function of a which is the control 
parameter in this setup. As expected, we find that E 3 
increases quickly when a decreases. For 1/4 < a < 1/2, 
E 3 is larger than E\ and E 2 by orders, thus we can safely 
ignore E\ and E 2 in this region. We get a one-particle 
system tunneling within an isolated double-well poten¬ 
tial, and should obtain two lowest energy states which 
are the bonding and the anti-bonding states separated 
by E 3 in energy. 



® [O 0 ] 9 [2n] 

FIG. 4. (Color online), (a) Four lowest energy levels as a 
function of <f>, with a — 0.45, E c — 0.05 Em, — 0.0001F m , 
and Jm — 5Em- (b) Wave functions of the two lowest energy 
levels with 4> = 0.50054>o; the dashed line is the potential 
function U~(0). 


V. FLUX QUBIT 

The analytic results shown in Fig. 3 are obtained for 
exactly half quanta of applied magnetic flux through the 
SQUID, where the double-well potentials are symmetric. 
In this case, the eigenstates are bonding and anti-bonding 
states for which the wave functions distribute equally at 
the two wells. If the applied magnetic flux deviate a little 
bit, the double-well becomes asymmetric. Then the low¬ 
est energy state should stay at the deeper well, while the 
other state should stay at the other well. We numerically 
solve the Hamiltonian in Eq. (j8|) to obtain the energy 
spectra. We apply both the finite difference method in 
real space and the plane wave expansion method in mo¬ 
mentum space, and find identical results. The four lowest 
energy levels of the spin-down particle as a function of 
applied magnetic flux 4> is shown in Fig. 4a. It can be 
seen that two lowest energy levels are isolated from other 
higher energy levels around <£> = 4>o/2. The two levels 
show opposite slopes as a function of <f>, indicating that 
the supercurrents carried by them have opposite circu¬ 
lating directions. To confirmation this, we calculate the 
wave functions of these two levels at <£> = 0.50054>o, as 
shown in Fig. 4b. As expected, we find that the two 
states stay at two different 0s, which account for con¬ 
tradicting super currents across the junctions. These nu¬ 
merical results demonstrate that the SQUID works effec¬ 
tively as a two-level system, with opposite supercurrents 
assigned to each level. We obtain identical energy levels 
for the spin-down particle. In fact, all levels are two¬ 
fold degenerate due to the spin degree of freedom. After 
checking the details of the wave functions, we find that 
the two lowest energy levels around <f> = 0.54>o can be 
described by an effective Hamiltonian, 

H= [A{$)t z +Bt x \v z , (17) 

where B « E%/2 determines the energy splitting at 
$ = $ 0 /2, and A($) » 8U ~^ |« = « o/2 ($ - $o/2) = 
[—7rJVl6<a 2 — l/4a 2 ](4>/4> 0 — 1/2). Here t X jZ are the 
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pauli matrices with the two supercurrent circulating di¬ 
rections as the basis states. The spin is conserved when 
the flux changes in the effective Hamiltonian Eq. (17), 
because the spin-flipping term is negligible comparing 
with the dominating energy scale B = E^/2. We obtain 
a typical effective Hamiltonian for a flux qubit, where the 
two basis states are represented by contradicting super¬ 
currents. 

We note that these numerical results are not sensitive 
to our present choice of junction parameters. In fact, 
similar results can be obtained for a range of parameters. 
From the quantum tunneling analysis, we find that the 
Josephson energy ratio a should be larger than 1/4 to 
keep the double-well potential, but smaller than 1/2 to 
achieve the isolation of one double-well. Our numerical 
results are consistent with the theoretical analysis, the 
two isolated levels described by Eq. (17) indeed exist 
within this range. 



FIG. 5. (Color online), (a) Four lowest energy levels as a 
function of <f>, with a = 1, E c = 0.05 Em, 5m — 0.015E m , 
and J m = 5 E m - The arrows show the spins of the two lowest 
energy levels, (b) Wave functions of the two lowest energy 
levels with <E> = 0.5014>o; the dashed line and the dotted line 
are the potential functions U~ (0) and U + (0), respectively. 


VI. COUPLING THE FLUX QUBIT TO THE 
MAJORANA QUBIT 

Looking back at Fig. 3, we notice that the spin-flipping 
hybridization energy E\ becomes important for large <a, 
since the other two hybridization energy become vanish¬ 
ingly small. Therefore, isolated low energy states should 
also exist. In this regime, however, the spin-flipping pro¬ 
cess plays the key role. Therefore, the spin is not con¬ 
served when changing the magnetic flux 4>. We confirm 
this picture by numerical simulations, and show the re¬ 
sults of the energy spectrum with a = 1 in Fig. 5a. We 
find two lowest energy levels which become energetically 
isolated from higher energy levels around <£> = 4>o/2. We 
further study the wave functions of these two levels at 

= 0.5014>o, and show the results in Fig. 5b. These two 
states exhibit contradicting supercurrent, therefore can 
be used as a flux qubit. Importantly, they have oppo¬ 
site spins, reflecting a spin hybridization as sketched in 
Fig. 5a. We check the details of the wave function of the 
states, and obtain an effective Hamiltonian for the two 
lowest energy levels around = 4>o/2, 

H = A(®)t z g z + B'r x a x , (18) 

where B' = E\/ 2. Unlike the pure flux qubit obtained 
in previous section, we arrive at a hybrid qubit where 
the flux qubit is coupled with the Majorana qubit (spin 
qubit). This hybrid qubit exists in the range of a > 0.7 
for the current choice of MBSs coupling 5 m . If 5 increases 
or decreases, the range of a should change accordingly. 

VII. DISCUSSIONS 

Our proposal is based on a topological nanowire 
Josephson junction. However, the idea could be naturally 
extended to other topological junctions where the 47r pe¬ 
riod Josephson effect exists. In particular, the quantum 


spin hall insulator in proximity to superconductivity is a 
promising candidate^. In addition, the 47 t Josephson pe¬ 
riod effect was also argued to be existing in topologically 
trivial Josephson junctions^, such as ballistic point con¬ 
tact junction. Our proposal should also work for these 
Josephson junctions. 

The proposed two-junction flux qubit could be real¬ 
ized with present experimental technologies. The key 
ingredient of the proposal, a topological superconduct¬ 
ing nanowire Josephson junction, has been reported in 
recent experiment s 32 ! 33 ! 40 . The growth of a conventional 
Josephson junction with controllable critical current was 
realized in several systems^. The superconducting cir¬ 
cuit which connects these two junctions into a SQUID 
could be achieved by the lithography technology, with 
controllable inductance and capacitance^. Moreover, 
these devices can be integrated easil y 1Q i 11 , therefore, the 
proposed flux qubit can be scaled to a large number. The 
typical value of the junction parameters are E m ~ 5/reU, 
E c ~ 0.1/ieV. The size of the SQUID should be in the 
range of micrometers with a vanishing self-conductance 
of L - lpH. The coupling energy of the MBSs depends 
on the length of the wire, with 5 ~ 1 neV for long wire 
and 5 ~ O.ljueV for short wire. The operation tempera¬ 
ture should be small to achieve quantum coherence, with 
an estimation of T ~ 50 mK. 

This two-j unction flux qubit can be measured by an 
external quantum circuit which detects the magnetic flux 
in the qubit accurately and sensitivel y 11 ! 12 . Meanwhile, 
it can be manipulated by applying an electro-magnetic 
pulse. When the frequency of the pulse fits the energy 
splitting between the two eigenstates of the flux qubit, 
the Rabi oscillation will drive a coherent rotation on the 
qubit stated. These operations on the two-j unction flux 
qubits do not require severer conditions than other types 
of flux qubits, thus should be experimentally realistic. 

We note that the operations on the flux qubit are 
based on electro-magnetic interactions. In this case, the 
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Major ana qubit should be untouched during the opera¬ 
tion since MBSs are charge neutral quasi-particles. This 
is correct for the pure flux qubit regime shown in Eq. 
(17). In this regime, the Majorana qubit is not coupled 
with the flux qubit, thus should remain the same under 
electro-magnetic perturbations. However, the situation is 
different for the hybrid qubit regime shown in Eq. (18). 
The coupling between the flux qubit and the Majorana 
qubit induces a strong correlation between flux state and 
the spin state, then the Majorana qubit can be manipu¬ 
lated by a rotation of the flux qubit. This manipulation 
method for the Majorana qubit might be useful for the 
topological quantum computation. 

VIII. CONCLUSION 

We propose a two-junction flux qubit using the 47r pe¬ 
riod Josephson effect. We reveal that this flux qubit 
shares all the merits of the previous two types of flux 
qubits, yet avoids their disadvantages. The proposed flux 
qubit is made of a small dc SQUID, formed by a con¬ 
ventional Josephson junction and a topological nanowire 
Josephson junction. The SQUID is built with a low in¬ 
ductance to reduce the external noises. In this regime, 
the system is fully described by a phase difference 0. The 
dynamics of the system resembles a spin one half particle 
moving in one dimension space, with 0 as the particle co¬ 
ordinate. The charging energy from the capacitance acts 


as the kinetic energy of the particle, and the Josephson 
energy acts as the potential energy. We investigate the 
Josephson energy of the system, and find that the 2tt 
period component in the nanowire junction can be an¬ 
nihilated by a destructive interference with the conven¬ 
tional junction. Then, the 47r period component in the 
nanowire dominates the system, bringing in a double-well 
shape potential. We study the low energy eigenstates of 
the system, and reveal that the double-well potential in¬ 
duces the isolation of two levels. These two levels have 
contradicting circulating super currents, thus can be used 
as a flux qubit. Furthermore, we find that a fine tuning 
of the SQUID will result in an coupling between the flux 
qubit and the Majorana qubit in the nanowire junction. 
This coupling might be important for quantum computa¬ 
tion. In summary, our proposal provides a two-junction 
flux qubit which has the merits of low inductance and 
large variability in junction parameters. 
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